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Abstract 



Confinement phenomenology characteristic of magnetically confined plas- 
mas emerges naturally from a simple sandpile algorithm when the parameter 
controlling redistribution scalelength is varied. Close analogues are found for 
enhanced confinement, edge pedestals, and edge localised modes (ELMs), and 
for the qualitative correlations between them. These results suggest that toka- 
mak observations of avalanching transport are deeply linked to the existence 
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of enhanced confinement and ELMs. 
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The introduction of the sandpile paradigm [|TJ- || into magnetized plasma physics (fusion 
|- [|12|], magnetospheric [HJ- |15], and accretion disk |16[- for recent reviews see Ref. 



fl9| , |20|| ) has opened new conceptual avenues. It provides a framework within which observa- 
tions of rapid nondiffusive nonlocal transport phenomena can be studied; recent examples 
include analyses of auroral energy deposition derived from global imaging [|15||, and of elec- 



tron temperature fluctuations in the DIII-D tokamak |L2[], both of which involve avalanching. 
Insofar as such phenomena resemble those in experimental sandpiles or mathematically ide- 
alized models thereof, they suggest that the confinement physics of macroscopic systems 
(plasma and other) may reflect unifying underlying principles. 

In this paper we present results suggesting that this unity may extend to some of the most 
distinctive features of toroidal magnetic plasma confinement: enhanced confinement regimes 
("H-modes"), edge localised modes ("ELMs"), steep edge gradients ("edge pedestal"), and 
their observed phenomenological and statistical correlations - for recent quantitative studies, 
see for example Refs. |2T|-p3[| and references therein. An important question is whether the 
L to H transition necessarily reflects a catastrophic bifurcation of confinement properties, or 
can be associated with a monotonic change in the character of the turbulence [^.Transi- 
tional behaviour resembling that observed in fusion plasmas has been found in other sandpile 
algorithms: for example in Ref [|J a local sheared flow region is imposed which reduces the 
characteristic avalanche length thereby affecting confinement; and in Ref |J, changes in the 
redistribution rule lead to changes in profile stiffness. We show that key elements of the ob- 
served phenomenology emerge naturally from a simple one-dimensional sandpile model, that 
of Chapman |p4}] , which incorporates other established models [|l|p5| as limiting cases. This 



centrally fueled (at cell n = 1) model's distinctive algorithmic feature relates to the local 
redistribution of sand at a cell (say at n — k) where the critical gradient z c is exceeded: the 
sandpile is instantaneously flattened behind the unstable cell over a length Lf, embracing 
the cells n = k — (Lf — 1), k — (Lf — 2), A;; and this sand is conservatively relocated to the 
cell at n = k + 1. Lf defines the length scale over which the most rapid redistribution occurs; 
in a plasma context this could be considered a proxy for turbulent correlation length or eddy 



size. In Ref. []24| the sandpile is explored for all regimes 1 < Lf < N for both constant and 
fluctuating critical gradient z c . Here we consider the dynamics of the more realistic case 
with random fluctuations in z c \ the system is robust in that once fluctuations are introduced 
in the critical gradient, the behavior is essentially insensitive to both their level and spectral 
properties [53], see also Ref. |TJ|. The limit Lf = 1 is the fixed point corresponding to the 



centrally fueled algorithm of Ref. |l] in one dimension. In the limit Lf = N (where N is the 
number of cells in the sandpile) the sandpile is flattened everywhere behind an unstable cell 
as in Refs. |25|JT9|| . A real space renormalization group approach |26f shows that the robust 



scale free dynamics for the limiting case Lf — N corresponds to a nontrivial (repulsive) 
fixed point (see e.g. Ref. ||). The essential result of Ref. ]2j| is that different regimes of 



avalanche statistics emerge, resembling a transition from regular to intermittent dynamics 
reminiscent of deterministic chaos. The control parameter is the normalized redistribution 
scalelength Lf/N which specifies whether the system is close to the nontrivial Lf = N fixed 
point. 

Height profiles for the sandpile with 512 cells, time averaged over many thousands of 
avalanches, are shown in Fig.l for three different values of the fluidization length Lf in the 
range 50 < Lf < 250. The sandpile profile shape, stored gravitational potential energy, and 
edge structure (smooth decline or pedestal) correlate with each other and with Lf. As Lf 
is reduced, the edge pedestal steepens and the time averaged stored energy rises; multiple 
"barriers" (regions of steep gradient) are visible in trace (a) and to some extent trace (b) of 
Fig.l. Time evolution of the sandpile for Lf = 50, 150, and 250 respectively is quantified in 
Figs. 2-4. The top traces show total stored energy; the middle traces show the position of the 
edge of the sandpile (the last occupied cell); and the bottom traces show the magnitude and 
occurrence times of mass loss events (hereafter MLEs) in which sand is lost from the system 
by being transferred beyond the 512th cell. Time is normalized to the mean inter- avalanche 
time At (proportional to the fueling rate). The sandpile is fueled only at the first cell, so 
that the great majority of avalanches terminate before reaching the 512th cell (these are 
classified as internal). While internal avalanches result in energy dissipation (recorded in 



the upper traces of Figs. 2-4), and may alter the position of the edge of the sandpile, they do 
not result in an MLE; there are corresponding periods of quiescence in the middle and lower 
traces of Figs. 2-4. Conversely the MLEs are associated with sudden inward movement of 
the sandpile edge, and in this important sense appear to be edge localised. However, MLEs 
and the associated inward edge movement are in fact the result of systemwide avalanches 
triggered at the sandpile center (cell n = 1). The character of the MLEs changes with 
Lf. In Fig. 2, where the mean and peak stored energy are greatest, the MLEs are similar 
to each other and occur with some regularity. The regularity of MLE occurrence in Fig. 3 
is less marked, the magnitude of the largest MLEs is greater than in Fig. 2, and there is 
greater spread in MLE size. This trend continues in Fig. 4, which also has the lowest stored 
energy. These effects correlate with the underlying dynamics of the sandpile. Figure 5 
plots the relation between average stored energy and Lf for the N = 512 system and much 
larger N = 4096 and 8192 systems (normalized to the system size N). The curves coincide, 
demonstrating invariance with respect to system size, with an inverse power law with slope 
close to —2 for Lf/N < 1/4, and a break at Lf/N ~ 1/4. These two regimes yield the 
quasi-regular and quasi-intermittent dynamics in Figs. 2-4 (see also the plot of avalanche 
length distribution against Lf in Fig. 8 of Ref. pi]). The parameter Lf/N is a measure 
of proximity of this high dimensional system to the Lf = N nontrivial fixed point. This 
determines both the apparent complexity of the timeseries in Figs. 2-4 and the underlying 
statistical simplicity described below, which is also invariant with respect to system size. 

There is systematic correlation between time averaged stored energy < E > and MLE 
frequency Jmle, as shown in Fig. 6. To obtain these curves, which are again normalized to 
system size, we have derived MLE frequencies using a standard algorithm previously used 
f22|j to assign frequencies to ELMs observed in tokamak plasmas in the Joint European Torus 
(JET). Since the sandpile often generates bursts of mass loss with structure down to the 
smallest timescales, which might not be resolvable under experimental conditions, we have 
followed Ref. [2~2"| in applying a (relatively narrow) measurement window of width 450Ar to 



obtain /mle- The correlation between < E > and fuhE is a noteworthy emergent property, 



furthermore Fig.6's characteristic curve is very similar to that of Fig.6 of Ref. [^TJ, which 
relates measured energy confinement to ELM frequency in JET. Energy confinement time 
r c can be defined for the sandpile by dividing the time averaged stored energy < E > by 
the time averaged energy dissipation rate < AE > (where AE is the energy dissipated in a 
single avalanche). The embedded plot of Fig. 6 shows r c against MLE frequency /mle- 

Finally, we explore the situation where there is a secular change in the redistribution 
algorithm: in Fig. 7, Lf decreases slowly, continuously, and linearly with time, from one 
constant value to another over a period encompassing many tens of thousands of avalanches. 
There is a corresponding time evolution of the energy confinement properties of the sandpile 
and of the character of the MLEs. Figure 7(top) shows total stored energy as a function of 
time as Lf changes from 250 at t = 4 x 10 4 to 50 at t = 1.15 x 10 5 , while ~ 10 5 avalanches 
occur: over a period of time corresponding to a few tens of MLEs, the system smoothly 
evolves from low to high confinement. This is accompanied by a gradual change in character 
of the time variation in the sandpile edge (position of last occupied cell, Fig.7(middle)) and of 
the MLEs (Fig.7(lower)), from large amplitude to small and from irregular to regular. Figure 
7 can perhaps be regarded as this sandpile's analogue of, for example, Fig.2 of Ref. [^2J or 
Fig. 2 of The essential point here is that the sandpile apparently freely explores phase 
space with changing control parameter Lf/N. Characteristic properties of the dynamics 
(whether quasi-regular or quasi-intermittent) and correspondingly, confinement properties 
(such as stored energy and MLE characteristics) smoothly follow changes in this parameter 
rather than exhibiting a sudden phase transition or catastrophe. 

By varying a single control parameter in the sandpile algorithm, we have shown corre- 
lations between: stored energy, confinement times, sandpile profile, sandpile edge structure, 
and the amplitude, frequency, and dynamical character of mass loss events. We have also 
seen how slow secular change in the control parameter produces a smooth evolution in con- 
finement properties. If a single control parameter analogous to Lf/N exists for tokamaks, 
it can in principle be found from experimental data by examining scaling with respect to 
system size as above. 
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The existence of such extensive tokamak-like phenomenology, emergent from a very sim- 
ple system, is a novel discovery. Insofar as the phenomeno logical resemblance is close, there 
is more to be learnt. A minimalist interpretation starts from the premise that the sandpile 
algorithm of Ref. |3J] provides a simple one-parameter model for studying generic nonlocal 
transport, conditioned by a critical gradient, in a macroscopic confinement system. Chang- 
ing the value of the single control parameter Lf then corresponds to altering the range in 
configuration space over which the transport process operates. It then follows from the 
results in the present paper that this may be the minimum requirement to generate those 
aspects of tokamak-like confinement phenomenology described. This is a significant conclu- 
sion, but one can consider a more far-reaching one. A possible maximalist interpretation 



attaches greater weight to recent observations p, 11,12] of avalanching transport in tokamaks 



and in largescale numerical simulations ||,|T0| thereof, and therefore regards the avalanch- 



ing transport that is built into sandpile algorithms as an additional point of contact with 
magnetically confined plasmas. One would then infer from the present results that toka- 
mak observations of avalanching transport are deeply linked to the existence of enhanced 
confinement and ELMs. 
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Captions 

FIG.l. Time averaged height profiles of the 512 cell sandpile for Lf = (a)50, (b)150 and 
(c)250. Inset: edge structure. 

FIG. 2. Time evolution of the 512 cell sandpile with Lf = 50: (top) stored energy, 
(middle) position of last occupied cell, (lower) magnitude and occurence of mass loss events. 

FIG.3. As Fig.2, for L f = 150. 

FIG.4. As Fig.2, for L f = 250. 

FIG. 5. Average stored energy versus Lf/N for sandpiles of N — 512, 4096, 8192. Energy 
is normalized to the Lf = 1 case (effectively to N 2 ). 

FIG. 6. Average stored energy versus MLE frequency, and (inset) r c versus MLE fre- 
quency for sandpiles of N = 512,4096,8192. Energy and MLE frequency are normalized as 
in Fig 6. 
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FIG. 7. Time evolution of (top) stored energy, (middle) sandpile edge position and (lower) 
MLEs, as Lf changes slowly and linearly from 250 to 50. 



10 



x 10 




520 



o 500 

Q_ 

CD 

■g> 480 

LU 

460 



6.5 



150 



n 



7.5 



8 



8.5 



9.5 



10 



^4 




x 10 



x 10 



LU 




6.5 





600 

CO 

o 

Q_ 

cu 400 

D) 
"D 
LU 



200 
6.5 

1500^ 



= 1000 
| 500 




6.5 



7.5 



8 



8.5 



7.5 



8 



8.5 



7.5 



8 8.5 
Time 



9.5 



9.5 



10 




x 10 



10 

4 




-2 -1.5 
log( L f / N ) 



-0.5 



500 



450 
400 
350 
300 

z 

^250 

LU 

V 

200 
150 
100 
50 



[ 
0. 



* 



+ 



* 



o 



o 



-©* 
o 

+ o 



+ 



* +o 

+ ' 



300 



200 



o 



100 



0L_ 
0.2 




0.4 0.5 

f MLE X N 



0.6 



*>-£H- 
i 



O* + * 



O ,* 



i ° + * i 



0.25 0.3 0.35 



0.4 

f MLE X N 



0.45 0.5 0.55 



0.6 



